Hyperfine structure and nuclear hyperpolarization observed in the bound exciton 

luminescence of Bi donors in natural Si 
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As the deepest group-V donor in Si, Bi has by far the largest hyperfine interaction and also a large 
1 — 9/2 nuclear spin. At zero field this splits the donor ground state into states having total spin 
5 and 4, which are fully resolved in the photoluminescence spectrum of Bi donor bound excitons. 
Under a magnetic field, the 60 expected allowed transitions cannot be individually resolved, but 
the effects of the nuclear spin distribution, —9/2 < I z < 9/2, are clearly observed. A strong 
hyperpolarization of the nuclear spin towards I z = —9/2 is observed to result from the nonresonant 
optical excitation. This is very similar to the recently reported optical hyperpolarization of P donors 
observed by EPR at higher magnetic fields. We introduce a new model to explain this effect, and 
predict that it may be very fast. 

PACS numbers: 78.55.Ap, 71.35.-y 
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Recent proposals [lHfJ to use the electron and nuclear 
spins of shallow donor impurities as qubits for Si-based 
quantum computing (QC) have led to renewed interest 
in these systems [5l4ll|. Most studies have focused on 
31 P, the most common donor in Si, with an / = 1/2 nu- 
clear spin. Most QC schemes involve enriched 28 Si, as 
this eliminates the 29 Si nuclear spin, but the removal of 
inhomogeneous isotope broadening [6| also enables an op- 
tical measurement of the donor electron and nuclear spin 
using the donor bound exciton (D°X) transition Q, and 
furthermore allows for the hyperpolarization of both spin 
systems at very low magnetic fields by resonant optical 
pumping 0. McCamey et al. have reported a dif- 
ferent effect in which P nuclear hyperpolarization can be 
achieved with nonresonant optical excitation in natural 
Si at high magnetic field and low temperature. 

Bismuth is the deepest group-V donor in Si, with a 
binding energy of 70.98 meV [121 ] ■ and is monoisotopic 
( 209 Bi), with a large 1 — 9/2 nuclear spin and a hyper- 
fine interaction of 1475.4 MHz, more than 12 times the 
117.53MHz value for 31 P [Hj]. While invoked in some 
QC proposals [4|], Bi has not been the subject of recent 
study. It is interesting to note that the Bi D°X in Si is 
described in the earliest studies of bound excitons (BE) 



in semiconductors 
tion since then [161 



14 



15j but has received little atten- 
This likely resulted from the scarcity 
of samples and, until now, their low quality. 



Recently 17] , Si:Bi samples have been grown from ul- 
trapure natural Si ( nat Si) using a floating-zone technique, 
for applications involving far- infrared lasers flij ]. Sam- 
ples from those same crystals are studied here, and show 
very reproducible D°X no-phonon (NP) photolumines- 
cence (PL) structure over a wide range of Bi concentra- 
tion. The spectra shown here are from a slice having a 
resistivity of 5.5f2-cm, mostly due to Bi, since the resid- 
ual B and P concentrations are estimated to be at least an 
order of magnitude less than the Bi concentration. The 



sample was mounted without strain in a high homogene- 
ity (0.01 %) split pair superconducting magnet dewar in 
Voigt configuration, with magnetic field B approximately 
along [100], and immersed either in liquid He or cold He 
gas. Above-gap excitation of 400 mW was provided by a 
1047 nm laser, and the collected PL was analyzed by a 
Bomem DA8 interferometer at a resolution of 1.8 /zeV full 
width at half maximum (FWHM). No effort was made to 
control either the excitation or the PL polarization. 

A level scheme indicating the expected PL Zeeman 
transitions between D°X and the donor ground state 
(D°) is shown in Fig [1] While this level scheme is sim- 
ilar to one introduced to explain the P D°X hyperfine 
structure in 28 Si |lfj,|ljj], it must be stressed that the P 
spectra were collected using PL excitation (PLE) spec- 
troscopy, in which the thermalizing initial states of the 
transitions were the D° states, whereas the present data 
is collected via PL spectroscopy, where the initial states 
are the D°X states, resulting in quite different thermal- 
ization behaviour. At zero magnetic field, the D° ground 
state is split into a doublet having total spins 5 and 4, 
and separated by 5 times the hyperfine interaction, or 
7377 MHz (30.51 /ieV) [l3j]. For nonzero B the D° states 
split according to the projection of the electron spin, S z , 
into two main branches, each of which has ten hyperfine 
subcomponents. The donor hyperfine levels are labelled 
from 1 1) to 1 20) in order of increasing energy (this order- 
ing only changes at extremelyhigh field) . As for all sub- 
stitutional donor D°X in Si [201] . in the ground state the 
two electrons must occupy the doubly-degenerate IsTi 
shell, forming a spin zero singlet, so the D°X splitting 
is determined only by the projection of the hole spin, 
J z (the nuclear Zeeman energy is small at these fields 
and is ignored, and the hyperfine coupling with the p- 
like hole is negligible). The six dipole-allowed D°X PL 
transitions are labelled from 1 to 6 in order of increasing 
energy, and each can be split into ten subcomponents by 



2 




B-Field (T) 



FIG. 1: (Color online) Zeeman level diagrams for Bi D° and 
Bi D°X indicating the origin of the PL structure from to 
2 T. J z labels the hole angular momentum projection in D°X, 
while S z and I z label electron and nuclear spin projections, 
respectively, in D°. The energy scale for D°X is compressed 
by a factor of 4. On the right the allowed transitions at 2 T 
are labelled from 1 to 6 in order of increasing energy, with 
each having 10 hyperfine subcomponents. 



the hyperfine interaction. 

The zero-field Si:Bi hyperfine doublet can be com- 
pletely resolved in all of our samples, as shown at the bot- 
tom of Fig. [51 The observed splitting is equal to the value 
expected from ESR [l3| , and the relative intensities are in 
good agreement with the degeneracies of the two states. 
The linewidths of 7.9 [ieV FWHM are slightly larger than 
the 5.7 /zeV width of the P D°X in nat Si. Given that this 
linewidth is also larger than the Bi hyperfine interaction, 
it is not surprising that individual hyperfine components 
cannot be resolved when a field is applied, although all 
60 hyperfine components are expected to be completely 
resolved in even a moderately enriched 28 Si:Bi sample. 
The hyperfine splittings are nevertheless evident in the 
Zeeman spectra even in nat Si, given the energy spread of 
the ten components spanning —9/2 < I z < 9/2, as can 
be seen in the three spectra in Fig. [5] taken at B = 2 T, 
the lowest field at which the six allowed PL transitions 
have no overlap. 

The top spectrum is at a temperature of T = 9K, as 
estimated from the width of the free exciton (FE) PL, 
and the six allowed transitions appear roughly rectangu- 
lar, as expected if the — 9/2 < I z < 9/2 states are equally 
populated (the total hyperfine energy difference between 
the —9/2 and 9/2 sublevels is much less than fc^T at 
any temperature used here). The observed spectra have 
been fit by calculating the transition energy of each of 
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FIG. 2: PL of the Si:Bi NP D°X transitions is shown at zero 
field, and at 2 T for 3 different temperatures. Vertical lines 
indicate relative strengths of individual hyperfine components 
obtained from the fit. The residual error of each fit is shown 
just below the spectrum. 



the sixty unresolved hyperfine components, and adding 
a line having a shape similar to the zero-field lineshape 
(but somewhat narrower) at that energy, with a relative 
intensity determined by D°X — > D° selection rules, ini- 
tial state thermalization, and the nuclear polarization. 
The known electron ^-factor is used in the fits, and the 
D°X hole ^-factors are adjustable parameters, together 
with the diamagnetic shifts [2l| and the temperature de- 
pendence of the band gap energy, all of which are op- 
timized across the spectra at various temperatures and 
B fields. An effective temperature, slightly higher than 
the nominal temperature, is obtained from matching the 
thermalization between the six main D°X transitions. 
The residual error of the fit is shown under each of the 
three spectra. 

A nuclear polarization term has to be included to ex- 
plain the skew of the observed spectra, especially at 
lower T, and given that the individual hyperfine com- 
ponents cannot be resolved, we make the simple as- 
sumption that the polarization per step of AI Z = 1, 
Pfit - [N(I Z + 1) - N(I Z )]/[N(I Z + 1) + N(I Z )}, is in- 
dependent of I z . The relative intensities and energies of 
the hyperfine subcomponents are indicated by the verti- 
cal lines under the PL components. It is clear that large 
nuclear polarizations are being produced, particularly at 
low T, where a substantial fraction of the nuclei are being 
polarized into the I z — —9/2 state. The nuclear polariza- 
tion results of the fits at B — 2 T, as well as at T = 1.5 K 
and B = 6 T are summarized in Table HI These nuclear 
hyperpolarizations under nonresonant optical excitation, 
in a direction opposite to what would be expected for 
the equilibrium nuclear polarization of the S z = —1/2 
branch, are very similar to the effect recently reported 
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TABLE I: Results of the PL fitting procedure for several B- 
fields and nominal T (Tbath.)- ^fit is the temperature deter- 
mined by the observed D°X thermalization, Pm is the nuclear 
polarization per step of AI Z — —1, P e is the expected equilib- 
rium polarization of D° electrons, and N(—9/2) is the frac- 
tion of Bi donors hyperpolarized into the nuclear spin state 
h = -9/2. 

B (T) T bath (K) T fit (K) P flt (%) P e (%) JV(-9/2)(%) 
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FIG. 3: PL spectra of the nat Si:Bi sample at low T and fields 
of 0, 3 and 6T. The excitation conditions and the intensity 
scale are the same for all three spectra. The inset shows a 
simple labelling scheme of the Bi and P hyperfine states for a 
discussion of the origin of the nuclear hyperpolarization. 

for nat Si:P at 8T M- 

McCamey et al. [9j explained this effect in terms of an 
Overhauser-like process, which we summarize with refer- 
ence to the inset in Fig. [3l The basis of their model is 
that the P electron spin thermalizes via the W processes 
at a temperature at or near the actual sample temper- 
ature T, whereas the R process equilibrates the |1) and 
1 3) states at some higher temperature T*. At low T and 
high B, where the equilibrium electron spin polarization 
is large, this will pump the D° population preferentially 
into |2). While this process may certainly play a role in 
the observed optical polarization, it is difficult to see how 
one could arrive at ab-initio estimates of the two tem- 
peratures, other than by using them to fit the observed 
polarization. 

We propose a different effect, inherent in the capture 
of FE to form D°X for all substitutional D° in Si, when 
B is high enough and T is low enough to generate signif- 



icant electron polarization. For substitutional donors in 
Si, the D°X ground state has two electrons with antipar- 
allel spins. At high field, both the D° electron and the 
FE electron will be well polarized into S z = —1/2, and 
the formation of D°X in its ground state requires the flip- 
ping of one of these electron spins, which remains ener- 
getically favourable since the D°X localization energy at 
the fields in question is still considerably larger than the 
electron Zeeman splitting. D°X having two S z = —1/2 
electrons might be an intermediate state in the formation 
of the ground state D°X, but this requires the captured 
electron to occupy the barely bound lsr 3i5 valley-orbit 
excited states, which lies ~4.3meV (~6.9meV) above 
the ground state for P (Bi) At T < 15 K these 

excited states are efficiently thermalized to the ground 
state within the D°X lifetime. In any case, formation 
of the ground state D°X requires an electron spin flip 
(AS Z = 1) which could occur via the spin-orbit effect. 

However, substitutional donors have another mecha- 
nism for achieving this electron spin flip during the cap- 
ture of polarized FE onto polarized D°. For P at high 
B and low T and without optical excitation, states |1) 
and |2) will each have nearly 50% of the total popula- 
tion (the energy difference between them being much less 
than ksT). Note that for P, states |2) and |4) are pure 
\S Z ,I Z ) states, whereas states |1) and |3), while tend- 
ing towards | — 1/2, 1/2} and 1 1/2 , —1/2), respectively, at 
high field, always have an admixture of the other com- 
ponent. D° in state |2) can only capture polarized FE 
via spin-orbit flipping of one of the electron spins. Those 
in state |1) can use the admixture of |l/2, —1/2) to flip 
the electron spin and form the D°X, while also flipping 
the nuclear spin for a total spin change of 0, and driving 
the nuclear population into I z = — 1/2, as is observed 0]. 
It is unlikely that the hole spin changes during the FE 
capture, but in any case the hole in the FE is already 
polarized into the J z — —3/2 state, so this would only 
increase the change in angular momentum. The same 
process will apply to Si:Bi, as shown on the left of the 
inset in Fig. [3j For Bi, only the 1 10) and 1 20) states are 
not mixed, and capture of polarized FE by D° in states 
|1) through 1 9) via this flip-flop process decreases I z by 
1, leading to a buildup of population in I z = —9/2. 

This model has testable consequences, namely, when 
the D° population builds up in the most favoured state 
(| 2) for Si:P and 1 10) for Si:Bi), capture of polarized FE 
to form D°X via the flip-flop process is no longer possi- 
ble. In Fig. [3] we show PL spectra of our nat Si:Bi sample 
on a wider energy scale at fields of 0, 3, and 6T, all at 
low T. The zero-field spectrum is dominated by the no- 
phonon line of the Bi D°X (BInp) and the transverse 
optical (TO) phonon replica of the boron acceptor BE 
(B^q). The optical phonon Raman line of the 1047 nm 
excitation laser is also observed. Even though the Bi 
concentration is much higher than that of B, the B PL is 
stronger, since the B BE has higher radiative quantum 
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efficiency than the Bi D°X which has a very short Auger 
lifetime [l6[- Note that B^ Q , the two-exciton B bound 
multiexciton complex (BMEC) is very weak at zero field, 
as expected for a sample containing ~ 10 15 cm -3 Bi, since 
the Bi capture most of the FE, keeping the FE concen- 
tration low. At 3T the absolute BInp intensity has de- 
creased, and the B^,q intensity increased, but even more 
noticeable is the large increase in the B 2 -^ intensity, and 
the appearance of Bip Q , indicating a large increase in the 
FE density (at zero field the lowest energy feature is the 
TO replica of the Bi D°X (Bi TO ), but at higher fields 
the Bi^ dominates). These changes are even more pro- 
nounced at 6 T. Between and 6 T the intensity of the Bi 
D°X PL relative to the B-related PL has decreased by a 
factor of 9. Almost exactly the same changes in relative 
intensity were observed between B and the P D°X in a 
nat Si sample containing only B and P. 

This supports our model for the origin of the opti- 
cally induced nuclear hyperpolarization, since once the Bi 
donors are hyperpolarized into the |10) state they can no 
longer capture polarized FE via this process, and hence 
the Bi PL decreases and the FE density increases, caus- 
ing the B BE and BMEC lines to increase in intensity 
(the electron(s) in the B BE and BMEC can all have 
S z = —1/2). Note that the lowest Zeeman state of the 
B BE and the D°X both have a single dipole allowed 
PL transition [20[, so the change in the relative PL in- 
tensities does not result from selection rules. The fact 
that B 1 has a single line at high B and low T while B 2 
has two and B 3 three is well understood [20]. The large 
increase in FE density made evident by the growth of 
the B BMEC lines at high field does not result from the 
fact that the FE are being polarized into a dark state 
with S z = —1/2 and J z = —3/2, since even at zero field 
the FE lifetime is completely dominated by capture onto 
donors and acceptors to form BE, a process much faster 
than the ms-long FE radiative recombination. 

We have been unable to observe transient polarization 
effects, leading us to suspect that the nuclear hyperpolar- 
ization might be too fast under the excitation conditions 
needed to observe PL for our low speed detector to follow. 
This surprising and potentially important possibility can 
be understood in terms of our model. Capture of a single 
polarized FE onto a D° in a mixed state should flip its 
nuclear spin with high probability, and under our exci- 
tation conditions this capture time is likely in the 10 to 
100 [is range. The ~100s polarization time observed by 
McCamey et al. @ could have resulted from much lower 
levels of optical excitation used in their experiment. 

In conclusion, we have resolved the zero-field hyper- 
fine splitting of the Bi donor in the PL spectrum of the 



NP line of the D°X transition in nat Si, and have ob- 
served a large nonresonant optical hyperpolarization of 
the nuclear spin under moderate to large B field, which 
is very similar to results recently reported for P At 
6T and 1.5 K, essentially all of the donors are pumped 
to the I z — —9/2 state. We propose a new model for 
the origin of this hyperpolarization, which is supported 
by changes in the overall PL spectra with magnetic field. 
Further experiments are needed to test the prediction 
that this hyperpolarization mechanism might in fact be 
extremely fast. PL in 28 Si:Bi would clearly resolve all 
60 D°X hyperfine components, and absorption or PLE 
measurements on the Bi D°X should result in a more 
direct measurement of the populations in all of the D° 
states. 
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